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ABSTRACT 

In order to expand the predictive capability of single- 
point turbulence closure models to account for the 
early-stage transition regime, a methodology for the 
formulation and calibration of model equations for the 
ensemble-averaged disturbance kinetic energy and en- 
ergy dissipation rate is presented. The calibration is 
based on homogeneous shear flow where disturbances 
can be described by rapid distortion theory (RI)T). The 
relationship between RDT and linear stability theory 
is exploited in order to obtain a closed set of mod- 
eled equations. The linear disturbance equations are 
solved directly so that the numerical simulation yields 
a database from which the closure coefficients in the 
ensemble-averaged disturbance equations can be deter- 
mined. 

INTRODUCTION 

Demands on the range of applicability of turbulence 
modeling are increasing, and with these increasing de- 
mands has come the need to develop models which are 
better sensitized to the transition process within the 
context of t he t raditional Reynolds-averaged turbulence 
modeling. Thus, a common mathematical framework 
linking the laminar regime, with its linear disturbances, 
and the fully turbulent regime, with its stochastic fluc- 
tuations, needs to be provided. 

As a first step tow F ard the development of such a 
mathematical framework. Thacker, Gatski and Grosch 
( 1999) studied the behavior of homogeneous, isotropic 
decaying disturbances. In turbulence modeling, this 
flow is used to calibrate the destruct ion-of-dissipation 
term in the modeled transport equation for the turbu- 
lent dissipation rate. In the linear disturbance case, 
the same functional form of the turbulence model was 
found but a new value of the model coefficient for the 


destrurtion-of-dissipation term was required. 

For consistency, disturbances are defined as devia- 
tions from the ensemble-mean in all flow regimes. The 
laminar regime is defined as the region of the How in 
which the ensemble-mean velocity is a stationary solu- 
tion of the Navier-Stokes equation. The disturbances 
in this regime are small enough in amplitude that their 
nonlinear' interactions can be neglected, and their evo- 
lution is completely predictable from their initial state. 
Traditionally, in laminar stability theory, disturbance 
fields are studied through the linear Orr-Sommerfeld 
equation, which describes the evolution of individual 
infinitesimal disturbance modes. Even when a quan- 
tity, such as the disturbance energy, is studied, it is the 
evolution of the instantaneous quantity rather than an 
ensemble average that is investigated. In contrast, the 
turbulent regime is defined as the region where the flow 
is subject to stochastic fluctuations, arising from non- 
linear int eractions, which render t he behavior of the dis- 
turbances unpredictable. In this case, the disturbance 
field is traditionally studied through (modeled) trans- 
port equations which describe the evolution of mean 
turbulent correlations. 

In this work, a set of (ensemble) mean disturbance 
transport equations, capable of describing the behavior 
of the lineai' disturbance fields, is developed. The ap- 
proach is founded on the observation that, even in the? 
laminar regime, every flow is subject to an inevitable 
uncertainty in initial conditions. Therefore, although 
each individual disturbance evolves deterministically, a 
probability distribution describing the initial (t = 0) 
energy distribution of the modes must be introduced 
for the calculation of ensemble-mean properties. This 
approach is similar to rapid distortion theory (RDT) 
in that it is based on linearized disturbance equations; 
however, the realm of application is different. RDT 
traditionally considers flows at higher Reynolds mini- 



hers in which the turbulence is fully developed and the 
effects of the molecular viscosity can be neglected, us- 
ing linearized equations to study the behavior of the 
disturbances under rapid (strong) distortion. The ap- 
proach taken here, on the other hand, considers linear 
disturbances in the early stages of transition where vis- 
cous effects must be taken into account- In addition, 
RDT is usually applied to short time evolution since 
in the turbulence case the nonlinear interactions cause 
sufficient growth of the fluctuations to render the linear 
approximation invalid after a few eddy turnover times. 
In this linearized disturbance case, no sucli limitation 
oil the time duration is encountered because viscous ef- 
fects result in the decay of the disturbance field at large 
times after a relatively small (initial) energetic growth. 

More recent extensions of RDT by Salhi, Cambon 
and Speziale (1997) have also exploited the connection 
with linear stability theory. They studied quadratic 
flows in a rotating frame to gain better insight into the 
dynamics so that a generalized stability criterion ap- 
plicable to turbulent flows could be developed. They 
also considered the effect on single-point, closure mod- 
eling - specifically the deficiencies in predicting elliptic 
flows. While the mathematical framework is similar in 
this study, the region of interest here is the early-stage 
transition regime. Nevertheless, this commonality fur- 
ther substantiates the basic assumption that a mathe- 
matical framework can be developed which will provide 
a set of transport equations capable of describing the 
flow (in a mean statistical sense) in the early-stage tran- 
, sition regime. 

In this study, the earlier analysis of Thacker et al. 
(1999) is applied to the case of mean homogeneous 
shear. Homogeneous shear flow is commonly used as 
a calibration flow for turbulence models because both 
turbulent transport and viscous terms can be removed 
from the transport equations for the turbulent correla- 
tions. The purpose of this study is to use the solution 
of the disturbance evolution equations for mean homo- 
geneous shear flow as a database in the calibration of 
the evolution equations for the ensemble-averaged dis- 
turbance correlations. As a first step, the focus is on a 
simple disturbance kinetic energy and disturbance dis- 
sipation rate (two-equation) description. In such a two- 
equation description of homogeneous shear, the terms 
in the kinetic energy equation are exact and require 
no modeling; whereas, in the dissipation rate equation 
both the production-of-dissipation and destruction-of- 
dissipation terms require modeling. The two closure 
coefficients associated with these terms are determined 
from the analysis presented here. Utilization of the dis- 
turbance evolution results as a reliable database is sup- 
ported bv the DNS results of Lee, Kim, and Moin (1990) 
who studied homogeneous shear flow at a high-shear 
rate. They showed that RDT results compared very 
well with the simulation results over the time period ex- 
amined. The results herein also show good agreement 
with the DNS results of Lee et al. (1990), and are found 
to apply at much later times due to the energetic decay 
of the disturbance field in the parameter range studied. 
Thus, this database will be used to provide insight into 
the asymptotic behavior of important dynamic vari- 
ables, as well as to provide the necessary information 
for the closure model calibration. The resulting closed 


disturbance dissipation rate equation can then be used 
in the formulation of a transition-sensitized turbulence 
model. 

LINEAR THEORY AND ENSEMBLE- 
AVERAGED CORRELATIONS 

In this section, transport equations for the ensemble- 
averaged linear dist urbance second- moments in homo- 
geneous shear flow are constructed. The theoretical de- 
velopment presented here parallels that of Townsend 
(1970) who described the structure of turbulence in a 
free shear flow as a product of the finite distortion of 
parcels of turbulent fluid. In this study, the relation- 
ship of RDT with linear stability theory (Speziale et al. 
1996, Salhi ct al. 1997) is expanded to include an analy- 
sis of t he transport equations for the ensemble- averaged 
disturbance kinetic energy and the disturbance energy 
dissipation rate. As in the turbulence case, such model 
equations require closure through the specification of 
closure constants. In the two-equation I\ — e formu- 
lation, only the disturbance dissipation rate equation 
contains modeled terms which have unknown closure 
constants. The homogeneous shear flow is used as a 
calibration flow for the production-of-dissipation and 
destruction-of-dissipation terms. 

In terms of dimensional coordinates (iq, i*2, i'3) in a 
fixed frame, the mean velocity is given by U\ — Sir 2, 
U> = b : 3 = 0, where the mean shear S = constant , 
and the disturbance velocity and pressure are denoted 
by Uj and ;3, respectively. For the problem of homo- 
geneous shear, it is convenient to work (Rogallo 1984) 
in a moving frame (fi , x 2 , * 3 ), in which the local mean 
velocity is zero and the coordinates are given by 

r 1 = i’l — ( St )i’2 

r 2 = x 2 

£3 = , (1) 

where t is time. The mean shear S and the viscos- 
ity v determine the time scale T = S~ i , length scale 
L = \/ uS ~ x , velocity scale U = \/ i/S , and pressure scale 
P = f 2 = 1 'S for the flow. This leads to the introduc- 
tion of dimensionless variables x 3 = x 3 /L (fixed frame), 
Xj = Xj/L (moving frame), t = t/T , u 3 = u 3 /U , and 
p = p/ P. The linear approximation is based on the 
assumption that the nonlinear terms in the dis- 

turbance momentum equations, can be neglected with 
respect to both the viscous term />V 2 u f and the mean 
shear term UjdU t /c)xj. This approximation remains 
valid as long as the amplitude of the disturbance ve- 
locity remains small compared with the velocity scale 
\/uS . 

An arbitrary solution of the incompressibility condi- 
tion and linearized Navier-Stokes equations in the mov- 
ing frame can be written as a linear superposition 

Ujlx'.t) = j d'k fj(k,t)e' k x ' 

p(x\f) = /d 3 kp(k,t)e' W . (2) 


:• the mode amplitudes /,(*./) 

anti p( k, t) satisfy 

fci/i +k' 2 h +h.h = 0 

<d> 

/, + /» + ik ip = -K 2 } , 

H) 

h + ik ’ 2 j, = -K 2 h 

(■■3) 

A. + ikxp = -K 2 h 

((>) 


with A> ~ k 2 - tk\ , A.'“ = Af 4* A^ 2 4 A 2 , and the wave 
numbers (Ai,A- 2, A?.) are constant in time. The solution 
for the velocity disturbance modes can be written in 
the form 

/,(k,f) = M u {k. t)/j(k.O)f . (7) 


is the probability distribution for disturbance second- 
moments, and the parameter ri, related to the variance 
of the distribution, will be given a physical interpreta- 
tion in what follows. This probability distribution is 
consistent with a A 1 low wave number behavior for the 
energy spectrum K{ k ) and, in the case of linear distur- 
bances in zero- mean-shear flow, leads to a decay 
law for the mean disturbance kinetic energy. 

The disturbance mode covariance is related to the 
energy spectrum tensor by 

{/, ( k, t)fj{\J)) = E (J (k, f )rf ?, (k 4 lb (Id) 

. so it follows from Eq.(7) that 


where tin* exponent 

r/( k. / ) = j K 2 dt 

= 1*2 r - A, h t 2 + {k,k,)t (8) 

arises from the viscous damping of all modes. The 
transfer matrix M has non- vanishing components 

A/n = Ah* = 1 


A 2 A 2 (0 — 0 O ) 

A/l2 = T-p 

k ' k H 

A 2 f(Ai — AaA'i 4 A 2 ) 

_ k;>k 2 ((-)o - 0) 

^ :y2 /4 

h H 

k [ A 3 t(k] f — A 2 A 2 4 k~ ) 

+ k 2 n A? 

wii h 

A/, = (fcf+fc?) 4 

00 = ,an_1 (fe) 

(-) = (an -1 ( 4 | (10) 

\kn ) 


E.jik.l) = A/,i(k, t) A! jm (k.t) 

xE 0 /„,(k) f - 2</lk ' ) , (14) 

where, as a result of (11). 

E a ,j{ k) = £ tJ (k,0) = (4 - A'oT’(k) (15) 

is the initial energy spectrum tensor. 

The energy spectrum tensor is the fundamental quan- 
tity from which all disturbance second- moments are de- 
rived by integration over wave number space. Of par- 
ticular interest are the disturbance stress tensor 

T,j(t) = j cl 3 k E,j(k,t) . (16) 

the disturbance kinetic energy K{t) = r„/2, and the 
disturbance dissipation rate 

s(<) = j (t k K 2 E„( k,f). (17) 

A physical interpretation can be given to the param- 
eter a, which appears in the probability distribution 
(12). by performing the integral in (17) at time zero to 
obtain £<> = 5/7 o/«. In terms of the dimensional distur- 
bance kinetic energy K = uSh\ and the dimensional 
dissipation rate i = uS 2 e the parameter a is given by 



The initial mode amplitudes /,( k,0) are elements 
of an ensemble with mean zero, and covariance 
(/ t (k, ())/;( 1, 0)). Assuming that the initial disturbance 
field is homogeneous, isotropic, and has finite energy 
A'o. the initial covariance is given bv the same expres- 
sion as that used by Thacker cl al. (1999) for the case 
of zero mean shear with no boundary 

</.(k.o)/,(l.O)) = <C(k + 1) (4 - ^-) 

x A'oT(k), ( 11 ) 

where 

**>- 1 ,i2) 


It will be shown that this quantity, a ratio of the dis- 
turbance time scale to the mean flow time scale, regu- 
lates the subsequent growth or decay of linear distur- 
bances in homogeneous shear flow. In the turbulence 
case, the variable >/ = SK/i also plays a critical role. 
Jongen and Gatski (1998) have shown that for turbu- 
lent homogeneous shear flow, r / reaches an equilibrium 
value which can be analytically connected to the vari- 
ation of the ratio of kinetic energy production to dissi- 
pation rate. 

MODELED DISTURBANCE TRANSPORT 
EQUATIONS 

The ensemble- mean kinetic energy and dissipation 
rate for lineal disturbances satisfy evolution equations 



which can he derived from the linearized Navier-Stokes 
equations (in the fixed frame). In this homogeneous 
flow, gradients of all ensemble-mean quantities vanish, 
and the (dimensionless) kinetic energy and dissipation 
rate equations reduce to 


K = -T 12 - £ = Vk - S, 

(19) 

i = V c -D c , 

(20) 


where Vk( — —Ty>) is the kinetic energy production, 


P = .,(/ <>». du, \ / dm du 2 \\ 

~\\i)xidx 2 / \dx, dx./J 

= -2 j d\[kik' 2 E..+IC 2 E l2 ] 

is the production-of-dissipation, and 

T>, = 

\dxjdxk dxjdx-k 

- ’ A "-"*' 


( 21 ) 


(22 


is the destruction-of-dissipation. From Eqs. (19) and 
(20), it is seen that only the disturbance dissipation rate 
equation requires closure models. In the standard A' — e 
formulation, the production- of-dissipat ion is modeled 
by 


V* = -T' £ 1~T12, 
A 

(23) 

which gives 


Cl --vJL-£L. 

ena £ Vk' 

(24) 

for the production-of-dissipation coefficient, 
struct ion- of- dissipation is modeled by 

The ile- 

Vc =c« 2 £, 
A 

(25) 

which yields 


c t2 = -v t 

(26) 


for the destruction-of-dissipation coefficient. 

Recall that the components E tJ of the energy spec- 
trum tensor are completely defined by Eqs. (7) through 
(15). With the solution of these equations, all the quan- 
tities of interest , including V t and V £ , can be obtained 
by carrying out the integration in Eqs. (16), (17), (21), 
and (22). These three-dimensional integrations were 
initially performed numerically using cubic and spher- 
ical grids. It was found that although the convergence 
was good at small times, at times greater than about 
105“ 1 , the results became strongly dependent on the 
grid resolution used (even for grids with as many as 
250 million points). The numerical evaluation of the 
integrals in Eqs. (16), (17), (21), and (22) is there- 
fore greatly simplified, and the accuracy significantly 


enhanced, by performing the radial integration analyt- 
ically. 

For the chosen probability distribution (12) the radial 
integrals in Eqs. (IG), (17). (21), and (22) are Gaussian 
and can easily be carried out analytically. Convergence 
tests show that the resulting angullar integrals had to 
be evaluated on grids with a spacing of 1/16 of a degree. 

In the next section the results of these computations 
are used to analyze the ensemble- averaged transport 
equations for the disturbance kinetic energy and dis- 
sipation rate. Of particular interest is the dissipation 
rate equation Eq. (20) which requires closure. It will 
be shown that closure models for V t and X> c , with con- 
stant closure coefficients, can be determined through an 
analysis of the long-time behavior of the disturbance 
correlations. 

RESULTS 

Figure 1 shows the evolution of the disturbance ki- 
netic energy for different initial values of ij o. For a 
small value of rfo (= 0.6), the energy decays monoton- 
ically in time. A transition between this monotonic 
decay and substantial growth of the disturbance energy 
occurs for values of rj o between 1.2, where the energy 
reaches a minimum and grows briefly, and 2.4, where 
the energy grows to a magnitude just equal to its initial 
value before decreasing. The temporal evolution shown 
in Fig. 1 for the moderate shear case 170 = 6, which 
corresponds to the equilibrium value for r/ in turbulent 
homogeneous shear flow, contrasts with the exponen- 
tial energetic growth (Tavoularis 1985) characteristic of 
the turbulent case. Here the disturbance kinetic energy 
grows to a little beyond three times its initial value. 
Therefore, if the initial disturbance kinetic energy is 
considerably less than the natural energy scale i/S , the 
linear approximation should remain valid for all times. 
Calculations have also been carried out for 170 = 12 and 
18. The evolution of A with time for these other 170 val- 
ues (not shown) is very similar to that for q 0 = 6. The 
only difference is that, as expected, the maximum value 
of K reached before the final decay is larger than for 
rTo = 6; for rj 0 = 12, K max = 6.95A'o, and for ?7o = 18, 
AVna.r = 10.82AV The initial value of t ] 0 = 18 is a high- 
shear case which is very close to the value used in the 
DNS study of Lee et al. (1990) (r/o % 17). Using the 
numerical simulation results, Lee et al. also compared 
with RDT and found very good agreement over the time 
interval examined (0 < t < 12). The high-shear limit- 
ing value 170 — » 00 corresponds to the case traditionally 
considered in RDT calculations (Rogers 1991), where 
viscosity is neglected. The absence of viscosity allows 
the energy to grow monotonically, rendering at some 
finite time the linear approximation invalid. 

Figure 2 shows qualitatively similar trends for the 
disturbance dissipation rate e(t) at the corresponding 
values of r/o = 0.6 and 6.0. For the small initial value 

= 0.G, the dissipation rate also decays monotoni- 
cally. At the moderate shear rate value 170 = 6, <r(t) 
grows after a brief initial period of decay, reaches a 
peak and then decays at a more rapid rate than the 
disturbance kinetic energy. Similar behavior was found 
for the higher shear cases of ijo = 12 and 18. 

Since the main focus here is on disturbances which 



can eventually grow (and lead to turbulence), tile re- 
sults for the initial value of 7/ 0 = 0.0 art 1 not of interest. 
In addition, while this weak shear case yields decaying 
disturbances, the validity of neglecting the nonlinear 
terms does come into question. As ( •amboii and Scott 
(1990) have pointed out, the contribution of the non- 
linear term, while small, may have a cumulative effect 
on the dynamics over a period of time. In addition, the 
linearized, weak shear case may be inconsistent because 
the product term containing the weak mean shear is re- 
tained while the nonlinear terms, possibly of the same 
magnitude, are omitted. These considerations dictate 
that subsequent, results will focus on the moderate- and 
high-shear cases. 

Fig. 3 shows the mapping of the anisotropy invariants 
//(= —b t j b,j /2 ) and / / / (= bikbkibtj /3), where b tJ (t) = 
Tij[t)f(2l\(t}) — S tJ /3 are the stress anisotropies. Re- 
call that initially, an isotropic distribution of the dis- 
turbance energy is assumed. The figure shows that af- 
ter t % 2, the (realizable) disturbance field migrates 
toward the two-component (2(9) state. After t % 10 
the disturbance field is 2(2 and monotonically evolves 
towards a one-component- (1G) state. Thus, the struc- 
ture of the disturbance field is significantly altered by 
t lie imposition of the mean shear. Also shown in Fig. 
3 is the invariant map trajectory of the DNS results of 
Lee et al. (1990). As can be seen, at early times (/ < 
2) both the DNS trajectory and that from the present 
calculations, with r/ o = 0 and 18. are in phase, and at 
later times, even though the DNS evolution lags behind 
the the current results, all of these evolve toward the 
same 1(2 state. While the results shown in Fig. 3 are 
for the r/o = 6 and 18 cases, the same trend toward a 1(2 
limit was found to hold for the case of 7/o = 0 . 6 . These 
disturbance field results are in sharp contrast to the 
equilibrium state reached for the turbulent case where 
the fixed point invariant values are III^ ~ 0.0043 and 
//oc » 0.064. 

Of particular importance for turbulence modeling are 
the closure coefficients C t \ and C s >. Figure 4 shows 
the evolution of these disturbance dissipation rate coef- 
ficients. Results for three initial values of rj o are shown 
which suggest that a limiting range of values at large 
times (and large 7 / 0 ) can be reached. For CV 1 there is 
a modest change in the value at t = 50 with changes 
in 7/0: for 7/0 = 0 . C.- t i = 2.23. for r/o = 12, (TVi — 2.08. 
and for r/ 0 = 18, C\ \ = 2.01. For (TV 2 , the values range 
about ±3*X from a mean of 2.46. Although it was not 
possible to determine analytically the limiting values 
for C c 1 and G £2 , numerical simulations for large values 
of r/o at large times showed that the limiting values were 
(\ 1 % 2.0 and (> t 2 ~ 2.5. These values are in contrast 
to the usual values obtained for turbulent closure mod- 
els where 6 £ i ss 1.45 and C\> lies in the range of 1.83 
to 1.92. 

In the turbulent, case, the destruction-of-dissipation 
coefficient was deduced from an analysis of the decay 
of isotropic turbulence (e.g. ( 'omte-Bellot and Corrsin 
1971). In the analysis of the decay of homogeneous, 
isotropic linear disturbances. Thacker et al. (1999) de- 
termined a value for CV 2 of 1.4 in a boundary-free case. 
This is in sharp contrast to the value of 2.5 found in 
this homogeneous shear flow. A 11 explanation for this 
difference may lie in the structure of the disturbance 


field. In the previous decay studies, the disturbance 
field was isotropic, and remained isotropic in the ab- 
sence of mean shear. However, as the results shown here 
indicate, the imposition of mean shear quickly produces 
ail anisotropic field ( see Fig. 3) whicli eventually drives 
the flow close to a 1(2 state (Fig. 3). Thus, a more 
relevant flow with which to form a comparison may be 
t he decay of anisotropic turbulence. Dakos and Gibson 
(1990) studied such a flow and deduced from the decay 
of the turbulence the destruction-of-dissipation rate co- 
efficient (\ 2 . Their anisotropic decay data yielded a 
value for ( ' c 2 of 2.18 which is relatively close to the 
value of 2.5 found here for (TV 2 Of course, in spite 1 
of fundamental differences in the dynamics associated 
with each study, it is interesting to note that the intro- 
duction of anisotropy in either disturbance field (linear 
or turbulent ) significantly increases the value of the co- 
efficient C t 2- 

CONCLUSIONS 

Closure models for the production-of-clissipation rate 
and destruction-of-dissipation rate in a linearized dis- 
turbance transport equation for the dissipation rate 
have been developed. The study lias extended the aj>- 
proach of Thacker et al. (1999) to mean homogeneous 
shear flow in which deterministic solutions of the lin- 
earized Navier- Stokes equation are combined with a 
probability distribution that accounts for the uncer- 
tainty in initial conditions, to obtain mean transport 
equations for disturbance correlations in the laminar 
regime. 

The temporal evolution of the disturbance kinetic en- 
ergy and the 1 disturbance dissipation rate wen* shown 
to depend on the magnitude of the initial value of the 
time scale ratio 7 /. The temporal evolution for both the 
kinetic energy and dissipation rate was monotonically 
decaying for values of r/o < 1.2. For larger values of 
r/o, both quantities initially decayed, then grew with 
time, and then subsequently decayed at large times. 
This is in contrast to the turbulent, homogeneous shear 
case where the quantities displayed exponential growth. 
Nevertheless, it was possible to find approximate val- 
ues for the dissipation rate coefficients C e i ^ 2.0 and 
C t 2 ~ 2.5, which were considerably larger than the cor- 
responding turbulent values. 
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Figure 2. Temporal evolution of disturbance dissipation 
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Figure 3. Invariant mapping of disturbance anisotropy 
tensor at various times for initial conditions: — ♦ — , 
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from DNS of Lee et al. (1990) at t = 2, 4, 6, 8. 10,12. 
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Figure 1. Temporal evolution of disturbance kinetic 

energy K for initial conditions: , rj 0 =0.6: , 
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Figure 4. Temporal evolution of closure coefficients C t \ 

and C £ 2 for initial conditions: , r/o =6; , 

rfo — 12; , t ] o = 18. 



